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We investigate the connection between Color Confinement and thermal Abelian monopoles popu- 
lating the deconfined phase of SU(2) Yang-Mills theory, by studying how the statistical properties of 
the monopole ensemble change as the confinement /deconfinement temperature is approached from 
above. In particular we study the distribution of monopole currents with multiple wrappings in 
the Euclidean time direction, corresponding to two or more particle permutations, and show that 
multiple wrappings increase as the deconfinement temperature is approached from above, in a way 
compatible with a condensation of such objects happening right at the deconfining transition. We 
also address the question of the thermal monopole mass, showing that different definitions give 
consistent results only around the transition, where the monopole mass goes down and becomes of 
the order of the critical temperature itself. 
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I. INTRODUCTION 



Color confinement is not yet fully understood in terms of the first principles of Quantum Chromodynamics (QCD). 
Models exist which relate confinement to the condensation of topological defects in the QCD ground state; one of 
them is based on dual superconductivity of the QCD vacuum 0, [2] ■ According to this model, color confinement is due 
to the spontaneous breaking of a magnetic simmetry, induced by the condensation of magnetically charged defects 
pH ' (e.g. magnetic monopoles), which yields a non- vanishing magnetically charged Higgs condensate. 
1— 1, The existence of a new phase of matter, in which quark and gluons are deconfined (Quark-Gluon Plasma), is a 
well defined prediction of lattice QCD simulations: the deconfined phase is under active experimental search in heavy 
ion experiments. However the physical properties expected for this phase are not yet clearly understood: the Quark- 
Gluon Plasma (QGP) is still strongly interacting above the deconfining temperature T c and its properties may be 
I more similar to those of an almost perfect liquid. 
— 1 , One of the hypotheses which have been put forward in the recent past is that QGP properties may be dominated 
by a magnetic component .3-5]. In Ref. Q| such magnetic component has been related to thermal Abelian monopoles 
. evaporating from the magnetic condensate which is believed to induce color confinement at low temperatures; moreover 
it has been proposed to detect such thermal monopoles in finite temperature lattice QCD simulations, by identifying 
them with monopole currents having a non-trivial wrapping in the Euclidean temporal direction 0, 0, 0j ■ First 
numerical investigations of these wrapping trajectories were performed in Ref. |6| and ML while a systematic study, 
regarding the deconfined phase of SU(2) Yang- Mills theory, has been performed in Ref. [8J. 

The definition of Abelian magnetic monopoles in non- Abelian gauge theories requires the identification of Abelian 
degrees of freedom: that is done usually by a procedure known as Abelian projection and relies on the choice of an 
adjoint field. Since no natural adjoint field exists in usual QCD, that implies some arbitrariness; a popular choice is 
to perform the projection in the so-called Maximal Abelian gauge (MAG). Results obtained in Ref. Q have shown 
that, as already well known for Abelian monopole currents in general, also the number and the locations of monopole 
currents with a non-trivial wrapping in the Euclidean time direction are quantities which depend on the choice of the 
Abelian projection. 

Despite that, the density and the spatial correlation functions of MAG thermal monopoles show a negligible de- 
pendence on the UV cut-off [|[ , as expected for a physical quantity. The temperature dependence of the monopole 
density, p, is not compatible with a (massive or massless) free particle behavior and is instead well described, in the 
whole range of temperatures explored, by a behavior p cx T 3 / (logT/ A e ff) 2 with A e // ~ 100 MeV, while the behavior 
p oc T 3 /(logT) 3 , predicted by dimensional reduction arguments, is compatible with data for T > 5 T c . This is in 
agreement with the picture of an electric dominated phase for Yang-Mills theories at very high temperatures, in which 
the magnetic component is strongly interacting Q . 

Moreover the study of density-density spatial correlation functions has verified the presence of a repulsive (attrac- 
tive) interaction for a monopole-monopole (monopole-antimonopole) pair, which at large distances is in agreement 
with a screened Coulomb potential and a screening length of the order of 0.1 fm. The above results have suggested a 
liquid-like behavior for the thermal monopole ensemble above T c Q and stimulated further research about the possible 
role of magnetic monopoles in the Quark-Gluon Plasma p^| - [l2j . 
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In the present paper, while aware of the problems related to the definition of thermal Abelian monopoles, we work 
on the hypothesis that those defined in the MAG may have a physical meaning and address a question regarding color 
confinement: if thermal monopoles in the deconfined phase are really related to the magnetic condensate responsible 
for confinement below the deconfincmcnt temperature, T < T c , is it possible to find clear signatures for their approach 
to condensation for temperatures slightly above T c ? Such approach would complement standard studies about the 
validity of the dual superconductor model, which look at the spontaneous breaking of a magnetic symmetry in the 
confined phase (i.e. the presence of a magnetic condensate) and at its restoration as T c is approached from below (]~3l — 

Since condensation is a phenomenon strictly related to quantum statistics, any signal of it should be searched 
for in properties of thermal Abelian monopole trajectories which reflect their nature of identical quantum objects. 
Such properties, as better explained in the following, are encoded in monopole trajectories which wrap two or more 
times in the Euclidean time direction and correspond to the permutation of two or more particles: no or very few 
multiple wrapping trajectories are expected if the system is very close to the Boltzmann approximation, while their 
number should increase as quantum effects becomes important, in a critical way close to the point where the particles 
condense. Such approach goes back to the seminal papers by Feynman (2lL [22J, where a path integral formulation 
was applied to describe condensation phenomena like the superfluid transition in 4 He. Since then, the path integral 
formulation of multi-particle systems has become an invaluable tool for analytical and numerical investigations of 
Bose condensation [23], [3l[ , and has been recently reconsidered in Ref . [24j , where it has been suggested to apply it 
also to the analysis of monopole condensation in QCD. 

A second question that we address regards the mass of thermal monopoles. Different strategies can be followed to 
determine some temperature dependent effective mass. For instance, one can look at the effective mass appearing in 
the description of multiply wrapping trajectories above T c or, in close analogy with the path integral formulation for 
non-relativistic quantum particles, one can obtain information about the mass by looking at the spatial fluctuations 
of wrapping trajectories. Of course it is not guaranteed that different strategies will lead to consistent determinations. 

The paper is organized as follows. In Section [XT] we recall a few technical details about the definition of thermal 
Abelian monopoles on the lattice. In Section HTT1 we study the distribution of multiply wrapping monopole trajectories 
in the deconfined phase of the SU(2) pure gauge theory and give an interpretation of them in terms of monopole 
condensation, determing a condensation temperature which is consistent with the known deconfinement temperature 
T c . In Section lIVI we address the question of the monopole mass. Finally, in Section [Vl we discuss our results and 
draw our conclusions. 



II. ABELIAN PROJECTION AND THERMAL MONOPOLES ON THE LATTICE 

The definition of Abelian monopoles requires the identification of a U(l) subgroup of the original gauge group. 
This is done by a procedure known as Abelian projection, which is assigned in terms of an adjoint field, i.e., in the 
particular case of SU(2) color gauge group, a vector field <f>(x) transforming as 

a-${x) -t G(x)(ff • $(x))G* (x) (1) 

under a local gauge transformation G(x) (a denote Pauli matrixes). If G fiv (x) is the field strength tensor of the 
original SU(2) gauge theory, then an Abelian tensor, known as 't Hooft tensor, is defined in terms of <fi(x) 

F^u = 4> ' G^ u — ~<j> ■ (D^(j) A D v <j>) 

= d ll {4>-A v )-dS-^)-- g h(9 tl 4>Ad^) (2) 

where <j>{x) = 4>(x)/\(j)(x)\. F^ v is an Abelian, gauge independent tensor, which in the gauge where (f>{x) = (0,0,1) 
takes the form 

= d^Af, — d u A^ . 

In that gauge, which is fixed up to a U(l) residual gauge freedom (cj> £ SO(3)/U(l)), the Abelian projection corre- 
sponds to taking the diagonal part of gauge links. 

In usual QCD there is no natural adjoint field, but several adjoint fields can be constructed in terms of gauge fields, 
typically a closed parallel transport (i.e. a path ordered product of gauge links). Alternatively, an implicit definition 
can be assigned by fixing = (0, 0, 1) and constant in one particular gauge: that defines (f>(x) in every other gauge. 



3 



An example is the so-called Maximal Abelian Gauge, where the gauge is fixed by maximizing the following functional 
with respect to gauge transformations: 

-Pmag = Rc tr [M n )^Ul(n) a a ] (3) 

//. n 

where Up(n) is a non- Abelian gauge link variable, i.e. an elementary parallel transport from lattice site n to n-\- jl. 
Fmag is proportional to the average squared diagonal part of the gauge links. 

In the gauge where the 't Hooft tensor is diagonal, Abelian link phases are extracted as follows: 

UJn) = u°Jd + iff-Up^ diag(e te " (n) , e-^ ( " } ) = V /' " = . (4) 

^K? + K) 2 

Abelian plaquettes, i.e. the lattice discretization of the 't Hooft field strength tensor, can then be constructed starting 
from the Abelian gauge link phases Op (n) : 

0^ = Op (n) + 6 V (n + A) - e„ (n + v) - 6 V (n) = dp6 v - d v dp . (5) 

Of course Abelian plaq uettes can be extracted also in a gauge where <fi(x) is not constant, even if through a more 
intricate procedure [25l|. obtaining equivalent results [25l [26] . as expected from the gauge invariance of 't Hooft tensor. 
Abelian monopole currents are identified by the standard De Grand- Toussaint construction |27| : 

rrifj, = —Suvpodjipa (6) 

Z7T 

where 8 ptT is the compactified part of the Abelian plaquette phase 

@p>i> = Op,u + 2-ku^ (7) 

and np V E N. 

Monopole currents form closed loops, since dprrin = 0. These loops may be either topologically trivial or wrapped 
around the lattice. Following the proposal of Ref. [4] , among the different currents we select those having a non-trivial 
positive (negative) wrapping around the Euclidean temporal direction, identifying them with thermal monopoles 
(antimonopoles) . 

This procedure identifies, at any given timeslice, the spatial positions of thermal monopoles with the points pierced 
by a current with non-trivial wrapping. There is some sort of ambiguity in this identification, since monopole 
currents extracted on the lattice are not strictly one-dimensional objects, but rather clusters of currents becoming 
approximately one-dimensional tubes in the high temperature phase, but staying anyway of finite thickness in lattice 
spacing units, because of ultraviolet (UV) noise which is present in the form of small monopole loops attached to 
the wrapping current. In our numerical setup we fix the spatial position of the monopole at the point where the 
wrapping monopole current is first detected by the current searching algorithm: thermal monopoles are thus spatially 
located with some random noise, which however has negligible systematic effects on the density and spatial correlation 
functions of thermal monopoles. 

Some of the monopole currents may wrap more than once, say k times, in the Euclidean time direction: in this case, 
as better explained in the following, we interpret them as sets of k monopoles (or antimonopoles) undergoing a cyclic 
permutation after a wrap in Euclidean time. The density and distribution of those fc-times wrapping trajectories is 
one of the subject of the present investigation. Also in this case, mainly due to the finite thickness of the wrapping 
current induced by UV fluctuations, ambiguities may be present for some specific configurations. One can imagine 
of situations in which it is hard to distinguish between two currents, each wrapping once, and a single current with 
a double wrap: that may happen, for instance, if the two currents overlap in some point and the answer may depend 
on the algorithm used to follow the current around the lattice. In our measurements we have chosen different current 
searching algorithms and we have taken the discrepancies between the different algorithms as an estimate of the 
systematic error linked to such ambiguities: such systematic error is always included in our determinations. 

In the following we shall present results obtained in the Maximal Abelian Gauge. It can be shown that, on stationary 
points of the MAG functional, the local operator 

X(n) = Y, [Up(n)<j 3 Ul(n) + U*(n - fi)a 3 Up(n - /})] (8) 

is diagonal. The maximization of the MAG functional for a given configuration has been achieved by an iterative 
combination of local maximization and overrelaxation (see Ref. [28|). stopping the algorithm when the average squared 
modulus of the non-diagonal part of X(n) was less than a given parameter ui [28|. We have chosen ui = 10~ 8 (see 
Ref. Q for more details). 
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FIG. 1: A possible contribution to the path integral representation of the partition function of 7 identical particles, in which 
the particles undergo a permutation made up of a 1-cycle, a 2-cycle and a 4-cycle. 



If we interpret the set of wrapping monopole trajectories, extracted from one gauge configuration of our Monte- 
Carlo sample, as one possible configuration of the Euclidean path integral representation of an ensemble of identical 
monopoles and antimonopoles at thermal equilibrium, then a trajectory which wraps k times before closing can be 
interpreted as a set of k monopoles (or antimonopoles) which permutate cyclically after going through the periodic 
Euclidean time direction. 

In the path integral describing N identical particles at thermal equilibrium, each possible configuration of the N 
particle paths contributing to the functional integral needs to be periodic, apart from a possible permutation of the N 
particles (the sign of the permutation is attached to the contribution if the particles are fermions) . That means that 
the configuration is not necessarily composed of N closed paths (that would correspond to the identical permutation), 
but is in general made up of M closed paths, with M < N: if the j-th path wraps kj times around the Euclidean 
time direction then X)f=i % = N anc ^ the configuration corresponds to a permutation made up of M cycles of sizes 
ki, fea, ... fejvf- I n Fig. [T] we report an example corresponding to a permutation of 7 particles partitioned into a 1-cycle, 
a 2-cycle and a 4-cycle. 

If effects related to quantum statistics are negligible, i.e. if the system is very close to the Boltzmann approximation, 
configurations corresponding to permutations different from the identical one are expected to have a negligible weight 
in the path integral, so that trajectories wrapping more than one time, corresponding to the exchange of two or more 
particles, are very rare. The number of trajectories wrapping more and more times is instead expected to increase 
as quantum effects become more important, and this should be especially true close to a transition associated with 
Bose-Einstein Condensation (BEC) (or with similar phenomena). 



In Table U we report the normalized average densities Pk/T 3 of trajectories wrapping k times as a function of 
temperature, determined for SU(2) pure gauge theory by use of the standard plaquette action and of the MAG 
Abelian projection. Data include both monopole and antimonopole trajectories and have been obtained by extracting 
monopole currents from samples consisting of O(10 3 ) independent gauge configurations for each value of T. The 
superscript a or b above the temperature value refers to two different values of the inverse gauge coupling, /3 = 2.70 



III. MONOPOLE CONDENSATION 



A. 



Numerical simulations 
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TABLE I: Normalized average densities pk/T 3 of trajectories wrapping fc times, as a function of T/T c for SU(2) pure gauge 
theory. The superscript a or 6 above the temperature value refers to two different values of the lattice spacing, 0.047 fm and 
0.063 fm respectively. Blank spaces indicate that no monopole current with the given wrapping has been observed in our 
statistical sample. 

and j3 = 2.60, corresponding respectively to lattice spacings a{0) ~ 0.047 fm and a{0) ~ 0.063 fm 1 . Simulations 
have been done at fixed spatial volume, while the temperature T — l/{N t a{P)) has been changed by varying the 
number of lattice sites N t in the temporal direction. In particular, at (3 = 2.70 we have made simulations on 64 3 x N t 
lattices, with Nt — 4, 5, . . . , 14, while at (3 — 2.60 we have chosen 48 3 x Nt lattices, with Nt — 4, 5, . . . , 10. Different 
spatial volumes have been considered in some cases for /3 = 2.60, in order to check for finite size effects. Confirming 
results reported for the overall monopole-antimonopole density in Ref. Q , data obtained for each p k show a reasonable 
independence from the value of the UV cutoff, even if small deviations are visible. 

Numerical data for Pk/T 3 are also reported in Fig. [5] What is apparent from the figure is that the relative weight of 
trajectories wrapping more than one time increases rapidly as T approaches T c from above. The number of monopoles 
or antimonopoles which are found in non-trivial cycles (i.e. k > 1) is less than 0.1% for T > 2.5 T c , meaning that the 
system is essentially Boltzmann-like at high temperatures. The same number goes to about 1% for T ~ 1.5T C and is 
well above 10% at the lowest temperature explored in our simulations. 

This is a clear qualitative signature for quantum statistics effects becoming more and more relevant as the critical 
temperature is approached, as expected if monopoles condense at T c . In the following we shall try to be more 
quantitative, but in order to do that we need to rely on some specific model. We shall first consider the analysis of 
a system of non-relativistic non-interacting bosons, following what performed in Refs. [22], [23[ as a starting point to 
describe the transition to superfluid helium. We shall then discuss the deviations from this very simple model which 
are expected for the monopole-antimonopole ensemble in QCD. 

B. Non-interacting Model for BEC 

We consider a system of N identical bosonic particles and their partition function Z = Tr(e~^/ T ). The trace 
can be taken over position eigenstates \x\ . . . a; at), but the correct states to consider in case of identical particles are 



The physical scale has been determined according to a(/3)A^ = R(/3)\(j3), where R is the two-loop perturbative /3-function, while 
A is a non-perturbative correction factor computed and reported in Ref. [29l |. We have assumed the values T c /A^ = 21.45(14) 29], 
Tc/Va = 0.69(2) [U and y/a ~ 430 MeV. 
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FIG. 2: Normalized densities pk/T 3 as a function of T/T c . 



proportional 2 to 

^2\x Pl ...x Pn ) 
p 

where the sum is over all the possible permutations P of the N particles. The partition function can be written as 

Z =^Y,J d 3 xi...J d 3 x N {x Pl ...x PN \e-P H \x 1 ...x N ). (9) 

We keep only the kinetic term K = p 2 /(2m) of the Hamiltionian H, discarding interactions and relativistic effects: 
we shall discuss these approximations in Section MI CI 

The assumption of absence of particle-particle interactions implies that the matrix element in Eq. ((9]) can be 
conveniently factorized according to the decomposition of the permutation into disjoint cycles. To clarify this point, 
consider the following explicit case involving 5 particles and a permutation composed of a 3-cycle plus a 2-cycle: 

(x 3 ,x 1 ,x 2 ,x 5 ,x 4 \e~ H/T \xi , x 2 , x 3 , x 4 , x 5 ) = 

(x 3:Xl ,x 2 \e-^ + ^ + ^/^ \x 1: x 2 ,x 3 )(x 5 ,x 4 \e-rt + rt/( 2mT > \x 4 , x 5 ) . (10) 

Also the integration in Eq. © can be carried on independently over groups of x variables belonging to the same 
cycle, so that each summand permutation can be factorized into a product of different contributions and the partition 
function can be rewritten as follows: 



^t£IR* en) 



Nl 

P k 



2 The normalization factor is 1/ y/N\ if the N coordinates are all different and changes otherwise. 
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where 

' d 3 Vl ... I d 3 Vk ( 2 / 2 , 2 / 3 ...2/ fc ,yi|e-^ + - + ^ )/(2mT) |yi,2/2...y fc ) (12) 



Zk 

is the contribution coming from one £:-cycle and n k is the number of different fc-cycles appearing in the cycle decom- 
position of the permutation P. It is straightforward to see that z k can be rewritten as 

z k {T) = J d 3 yi ...J d s y k { Vl \e-^\y k ) ...{y li \e~^r\y 2 ){y 2 \e-^r\y 1 ) 

= [ d 3 y 1 (y 1 \e- k ^\y 1 )=z 1 (T/k) (13) 



where zi(T/k) is the partition function for a single particle at a temperature k times lower. Notice that the whole 
derivation would be unchanged if the particles move in an external potential, but still in absence of particle-particle 
interactions. 

The partition function for a non-relativistic free particle in a box is well known: 

Zl {T)=V/\ 3 (14) 



where V is the available volume and A = y^^n/ (mT) is the De Broglie thermal wavelength. There is an approximation 
involved in the result above: the spacing between the quantum energy levels of the particle must be negligible with 
respect to T, so that the sum over energy levels can be replaced by an integral. That condition can be written 
equivalently, for a cubic box, as A <C V 1 ^ 3 , i.e. the thermal wavelength must be neglibible with respect to the size of 
the box; a similar condition applies if the particle is constrained in a three-dimensional periodic torus, as happens in 
our lattice setup. For a given volume, this condition is violated at low enough temperatures. 
Combining Eqs. ([13")) and (IT4l we obtain 

Zk{T ) = z 1 (T/k) = ^ m (15) 

but the condition for its validity is now \fk A <C V 1 ! 3 . Such condition can be violated, for given V and T, if k is 
large enough. To better understand this condition, think of the loop constructed with the positions of the k particles 
belonging to the cycle; the integral in Eq. (|12l) can be rewritten in terms of integration over the relative coordinates 
of neighbooring particles plus an integration over the center of mass of the loop: only if the size of the loop, which 
is of order \fk A, is negligible with respect to the size of the box, the integration over the center of mass can be 
approximated by V. 

In the following we shall substitute Eq. (|15p into Eq. (|11[) . assuming that the contribution from permutations 
having very large cycles is suppressed: that coincides with the approximation made in Ref. [23j. The approximation 
is expected to work well, apart from very close to the condensation temperature. Data reported in TableUgive us the 
opportunity to check directly if this assumption is reasonable for our ensemble of monopole trajectories. We see that 
indeed large values of A; are exponentially suppressed, so that only cycles with small k can be found, apart from very 
close to T c . Therefore we shall assume the validity of Eq. (fl"5|) . being aware that violations are expected close to T c : 
such violations are part of the finite size effects arising close to a critical point. 

With the assumption above, the partition function in Eq. (Ill|) becomes 

s-siEnGW- s ^-■*-»»>n(£ i kr (16) 

P k v ' ni,fls...Tliv k v ' 

where we made explicit the partition into cycles (n\ 1-cycles, n 2 2-cycles, . . . n k fc-cycles) of each permutation P. The 
factor C(ni, n 2 . . . n/v) accounts for the multiplicity of each partition {m, n 2 . . . n^}: 

N\ 

C{n 1 ,n 2 ...n N ) = - — - , (17) 

l ni 2™ 2 . . . n\\n 2 \ . . . 

as it follows from a simple calculation. 

The exact computation of the sum in (fTS]) is not simple, due to the presence of the constraint kn k = N 
on the possible partitions, which must sum up to the total number of particles. This difficulty can be avoided by 
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relaxing the constraint, i.e. by switching to the Grand Canonical formalism. The Grand Canonical partition function, 
Z = Tr(e~( H ~ ,J ' N " T ), can be easily computed starting from Eqs. (fTBf and (fTT|) : 

*~ 2^ i"i2"2...„ 1 ! Tl2 !...lll A 3 P/ 2 J ~H P lA 3 fc 5/2 7 

ni,?i2... fe fe 

and from the Grand Canonical partition function the average density of fc-cycles easily follows: 

_ (n k ) e-» k 

where we have defined the dimensionless chemical potential p = —p/T, which satisfies the constraint p > (i.e. 
M<0). 

The expression for the total particle density is then 

N A i ~ e -^ fc 

y- p = ^ kpk = #T,W*- (20) 

k=l k=l 

Of course last expression must be equal to the usual result for the density of an ideal non-relativistic boson gas which 
is obtained when working in momentum space: that can be easily verified by recalling one of the integral expressions 
for the polylogarithm Li s (z) = z k /k s , in our case in particular: 



E 



dx-r^ ■ (21) 



fc 3 / 2 Jtt J n ei 1 e x - 1 

k v u 

The expression for the total density is bounded by an upper limit which is reached for p = 

1 ^ i 2.612 

fe=l 

For larger densities or, at fixed density, for lower temperatures (larger A) our description is inadequate and it is 
necessary to allow for macroscopic cycles or equivalently, when working in momentum space, for a macroscopic 
number of particles occupying the ground (zero momentum) state: that corresponds to Bose-Einstein condensation. 
The critical temperature (density) is signalled by the vanishing of the chemical potential. 

If we had a set of trajectories in configuration space sampled from the path integral representation of an ideal 
non-relativistic Bose-Einstein gas at various different temperatures, we could measure the distribution of fc-cycles, i.e. 
the densities pk, and then, by fitting the expected dependence given in Eq. (fT9|) . we could determine the chemical 
potential fi as a function of T, in order to find numerically the critical temperature Tbec a t which /t vanishes and 
Bose-Einstein condensation begins. 

The advantage of treating the condensation of the ideal gas in configuration space rather than in momentum space 
should be obvious at this stage: we have a direct connection with the information about the monopole ensemble 
which can be retrieved by lattice QCD simulations. The aim of our analysis is indeed to try a similar thing with our 
ensemble of monopole trajectories, i.e. to determine if some critical temperature Tbec exists at which the analogous 
of the chemical potential vanishes, and what is its relation with the deconfinement temperature T c . There are of 
course various caveats related to particle-particle interactions, that we shall discuss in the next subsection. 



C. BEC of thermal monopoles 



The ensemble of thermal Abelian monopoles that we are investigating is surely different from an ideal gas of non- 
relativistic identical particles. The analysis of Ref. 0] indeed has shown the presence of particle-particle interactions, 
which are attractive in the monopole-antimonopole case and repulsive in the monopole-monopole case. The presence 
of those interactions interferes with many steps of the above derivation for the densities pk, for instance it is not 
possible to rewrite the contribution of a single fc-cycle as in Eq. (fT5)) . 

One of the effects of interactions is visible in the particle density around a given thermal monopole, which is 
suppressed (enhanced) for monopoles (antimonopoles) with respect to the average density. In the monopole-monopole 
case the suppression extends, close to T c , over distances of the order of 1 fm: as we shall discuss in the next Section, 
such distance is larger than the typical thermal wavelength A which can be measured on trajectories wrapping only 
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FIG. 3: Fit of the densities pk according to e /k 5 ^ 2 (dashed line) and according to e ^ k (solid line) for two values of the 
temperature. 
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TABLE II: Chemical potentials obtained for different temperatures by fitting the densities pk according to pk oc e ^ k /k a , for 
different values of a. 

one time. As a consequence, interactions should clearly affect monopole trajectories corresponding to fc-cycles with 
k > 1; the typical size of these loops should be \Vk in the non-interacting case, but is expected to be larger due to 
monopole-monopole repulsion: that corresponds indeed to the outcome of the analysis to be presented in the next 
Section. 

Apart from the interactions among particles belonging to the same loop, one should also take into account loop- 
loop interactions, which can be both repulsive or attractive (a monopole /c-cycle will repel (attract) a monopole 
(antimonopole) fc'-cycle), so that it is also not possible to factorize the contribution of each permutation into the 
contribution of different fc-cycles, as in Eq. (fTT]) . 

Even if we have some partial information about interactions from the analysis of Ref . [8] , taking them properly into 
account, in order to repeat the description of BEC along the same lines reported above for the free case, is not an 
easy task. 

Instead of looking for an exact solution, let us try to understand which similarities and differences should be 
expected for the dependence of the density pk on k, with respect to the free case reported in Eq. (TTI)|) . On general 
grounds one may expect some finite free energy cost needed to add one particle to a fc-cycle, i.e. to go from k to k+ 1, 
playing the role of an effective chemical potential, plus some interaction dependent contribution, so that: 

Pk = e-» k f(k) (23) 

where f(k) is some function decreasing less than exponentially with k. We shall refer to fx as chemical potential in 
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FIG. 4: Chemical potentials reported in Table HT1 for a — and a — 2.5 and two different lattice spacings, together with a fit 
of them according to Eq. (|24l) . 
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TABLE III: Results of the fit of chemical potentials according to jx = A (T — Tbec)" , for various different values of a. The 
lowest temperature, T = 1.017, has always been discarded from the fit. 

the following, however it should be clear that in the present context, in which we are dealing with a neutral plasma 
of monopoles and antimonopoles, it should be better regarded as a parameter for an effective description of the 
distribution in the number of fc-cycles. 

The result obtained in the free case, f(k) = (A 3 ^ 5 ^ 2 )" 1 , can be modified by interactions in various ways. In the case 
of a dilute hard sphere gas model, reported in Ref. [23|], f(k) is the same as in the free case plus corrections of order 
r/A, where r is the sphere radius. Part of the interactions could be also taken into account by an effective dynamical 
mass, as done by Feynman for the study of 4 He. In general we may expect a leading contribution /(fc) cx l/k a , 
where a could be different from 5/2 (this is also the way in which relativistic effects should show up). In any case, 
the approach to condensation should be signalled by the vanishing of the chemical potential, i.e. at the condensation 
point large fc-cycles should cease to be suppressed exponentially in fc. This is exactly what we want to check on 
our data reported in Table |TJ As we shall discuss soon, the outcome is that the vanishing of the chemical potential 
happens at a point which is compatible with T c and that this result is remarkably stable for various different choices 
of the function /(fc). 

We have tried to fit our data for the normalized densities Pk/T 3 according to Ae~^ k /k a . In order to obtain 
reasonable values for x 2 /d.o.f. we had to take into account only data with k > 2 — 3 for T < 1.2 T c and with fc > 1 
for higher temperatures. If the a parameter is left free, the x 2 is minimized for a around 2 for most of the explored 
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temperatures, however reasonable fits are obtained for a larger range of a values 3 . A few examples of such fits are 
shown in Fig. [31 while in Table [TT] we report results obtained for the chemical potential for various values of a going 
from to 3. 

As a second step, we have tried to fit the values obtained for the chemical potential according to a critical behavior: 

p = A (T - T BEC r' (24) 

which predicts p to vanish, as the condensation temperature Tbec is approached from above, with a critical exponent 
v' . Such functional dependence describes reasonably well data reported in Tabic [Til f° r each different value of a, if 
the lowest temperature, T = 1.017 T c , is discarded: this is expected, since the vanishing of p implies the appearance 
of fc-cycies with arbitrarily large k, i.e. cycles of permutating particles with arbitrarily large spatial extension, 
corresponding to a diverging correlation length; this is not possible on a finite lattice hence the dependence Eq. (|24|) 
must be valid not too close to the condensation temperature. One might wonder whether such finite size effects are 
also present, on the lattices used for our investigation, for higher values of T: this issue is discussed in detail in 
Section llll Dl where we show that this is not the case already for T — 1.052 T c . 

We report in Table HTT1 the results for Tbec an d v' obtained from fits limited to a temperature range 1.052 < 
T/T c < 1.582; values do not change within errors if the range is slightly changed (but not including the lowest 
temperature T = 1.017 T c ). A few examples of such fits are also reported in Fig. [4] The values obtained for x 2 /d.o.f. 
are reasonable, taking into account that we have mixed data obtained for different values of the lattice spacing and 
that scaling violations, even if small, may be present; the ansatz given in Eq. (|24p seems to work slightly better for 
the chemical potentials obtained in the case a = 2.5, i.e. the value expected for free bosons. 

As can be inferred from Table HTT| the dependence of Tbec and v' on the exponent a, used to deduce the chemical 
potentials from the densities pk, is very mild or even negligible within errors. That makes the main result of our 
analysis stronger and less dependent on the way we deal with interactions: Tbec is compatible within errors with 
the critical temperature T c at which the SU(2) pure gauge theory deconfines. Therefore we have proved that, as we 
proceed from higher to lower temperatures, the onset of confinement is associated with the condensation of thermal 
monopoles present in the deconfined phase, as expected in the dual superconductor scenario for color confinement. 

If monopole condensation indeed coincides with deconfinement, a second important question regards the relation 
of the critical exponent v' with the critical exponents of the second order SU(2) deconfinement transition, which 
belongs to the 3d Ising universality class. As clearly visible from Table Unl the value obtained for v' is compatible or 
slightly larger than the critical exponent related to the diverging behavior of the correlation length at the transition, 
£ ~ (T — T c ) _!y , which for Ising 3d is v ~ 0.63: in the following we shall try to give a possible interpretation of this 
fact. 

The increase in the correlation length £, observed when approaching the critical condensation temperature from 
above, may be related to the appearance of larger and larger fc-cycles, corresponding to groups of permutating particles 
which extend over larger and larger spatial extensions. We may assume that £ grows proportionally to the typical 
spatial extension of fc-cycles. For a distribution of fc-cycles pk ~ e~^ fc /fc a , the average value of fc is (fc) oc 1 /p; therefore, 
if we knew how the typical spatial extension of a fc-cycle depends on fc, the dependence of £ on /t would be easily 
obtained. In the free boson case, the typical spatial extension of a fc-cycle is oc y/k, corresponding to a random walk 
behavior. However the size is expected to grow faster with k because of particle-particle repulsive interactions: as a 
limiting case, opposite to the random walk behavior, the tyical size of a fc-cycle may grow linearly with fc if interactions 
induce ordered linear structures in the set of permutating particles (e.g. the particles are typically disposed on some 
closed smooth line). In general we expect £ cx 1/p^, where we have introduced the exponent u> which dictates how 
the size of fc-cycles grows and is expected to be in the range 1/2 < to < 1: lower values of ui could be possible if the 
permutating particles are more closely packed than in the non-interacting case, but this is not expected in presence 
of repulsive interactions. Since £ ~ (T — T c )~ v , we expect p ~ (T — T c ) u /^ , hence v' — v/u>. 

From Tabic Unl we conclude that v' is compatible with v, as expected if oj ~ 1: this suggests that the typical sets 
of permutating thermal monopoles contributing to the path integral form linearly ordered structures, which are likely 
induced by interactions. Therefore the typical cluster of monopole currents corresponding to a fc-cycle contribution 
should lie on a sort of time oriented closed surface (see Ref. for a recent discussion about the possible appearance 
of similar structure right above T c ) . 



3 Notice that a value a > 2 is to be preferred for the asymptotic large k behavior, since in this case the total number of particles ^ kp^ 
is convergent even for fi = 0, thus claiming for a macroscopic number of particles occupying the ground state. 
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TABLE IV: Normalized densities for T — 1.052 T c (f3 = 2.60 and N t = 10) and various spatial lattices L s . The chemical 
potentials obtained for a — 2.5 are also reported. For the explorable number of wrappings, finite size effects are negligible, 
within errors, for L s > 48. 
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FIG. 5: Dependence of the normalized densities pk/T 3 on the spatial lattice size for /3 = 2.60 and Nt = 10 (T/T c = 1.052). 



D. Finite size effects 



As discussed above, it is natural to expect that the density of trajectories with a high wrapping number may 
be affected by the finiteness of the system in the spatial directions. That may be a source of systematic error in 
our analysis, in particular for the temperatures closest to the transition, where the density of multiple wrapping 
trajectories is higher. In fact, if the typical cycle extension is really linked to the correlation length of the system, 
such effects are another way of looking at the saturation of £ as the critical temperature is approached on a finite 
lattice. Our data have been obtained on two different combinations of lattice size and bare coupling, i.e. L s — 48 at 
/3 = 2.6 and L s — 64 at /3 = 2.7: in both cases the physical spatial size is approximately the same, T c L s a((3) ~ 4.5. 

Therefore, in order to check for finite size effects, we have considered the lowest temperature included in our 
analysis, i.e. T = 1.052 T c (corresponding in our numerical setup to /3 = 2.6 and Nt = 10), where we have repeated 
the determination of the densities of multiple wrapping trajectories and of the chemical potential fi on various different 
spatial lattice sizes, L s = 32, 40, 48, 56, 64, corresponding to T c L s a((3) ranging from 3 to 6. 

Results obtained for the densities of multiple wrapping trajectories are reported in Table IIVI and in Fig. [5] In 
Table HVl we also report the determination of the chemical potential jl in the free boson approximation {a — 2.5) on 
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TABLE V: Monopole masses, determined for different temperatures according to the definitions given in Eqs. (|25[) and ((28 
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FIG. 6: Monopole masses as a function of T/T c , as reported in Table IVl 



the different lattice sizes. 

Finite size effects are clearly visible on the smallest lattice and result in a stronger suppression of trajectories with 
a large number of wrappings, or equivalently in a higher value of fi. However results are compatible within errors for 
the three largest lattices, L s = 48, 56, 64. We take that as evidence that finite size systematic effects are negligible, 
within our current statistical uncertainties, for spatial lattice sizes aL s > 4.5 Tq 1 and for temperatures T > 1.052 T c . 
Such systematic effects therefore do not affect our analysis presented above. 



IV. MONOPOLE MASS 



We have followed different strategies in order to determine a temperature dependent effective monopole mass. As 
we shall discuss, different definitions lead to different results, so that no unambiguous determination can be achieved; 
nevertheless one can identify some common features. 

The first strategy is to consider the fit of the densities of wrapping trajectories according to the non-relativistic 
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TABLE VI: Average squared fluctuations measured as reported in Eq. (|27[) . as a function of T/T c , for trajectories having 
wrapping numbers k — 1,2, 3. 



free boson prediction, Eq. (TT"9"f , and extract from the coefficient 1/ A 3 an estimate for an effective dynamical mass 



2tt 

m BEC = yTTa ■ ( 25 ) 



If we consider the same fits used to compute the chemical potentials reported in Table ITT1 for a = 2.5, we obtain the 
masses reported in the second column of Table [V] and also in Fig. [B] The dependence on the lattice spacing seems 
negligible within errors. If we neglect the determination obtained for T = 1.017 T c , where as discussed in previous 
Section finite size effects could be significant, we see a clear trend for an effective monopole mass rapidly decreasing 
from values of the order of 20 T c for T ~ 1.5 T c down to tobec ~ T c at T c . Of course when meEC is of the order 
of T c the non-relativistic approximation breaks down, however even a linear fit taking into account only values for 
T > 1.095 T c , where wbec is already a few times T c and relativitic corrections can be neglected, leads to predict 
hibec ~ T c at the transition. 

A different strategy is to relate the monopole mass to the spatial fluctuations of monopole trajectories: those should 
be more and more damped as the mass increases. Consider again the path integral of a non-relativistic particle of 
mass m and the average squared spatial fluctuation of the particle periodic path, defined as 

Ax 2 = T f ' dt((x(t) - x(0)) 2 ) , (26) 



where the average is taken over the path integral distribution. That is directly related to the mass of the particle. 
A simple calculation shows that for a free non-relativistic particle one has Ax 2 = 1/ (2mT) , while for a particle in a 
harmonic potential with elastic constant muS 2 one obtains: Ax 2 = -6/(/3mw 2 ) + 3coth(w/3/2)/(mw). 
On the lattice the mean squared monopole fluctuation is defined as 

a- 2 A* 2 = ±X> 2 , (27) 
t=i 

where di is the squared spatial distance (accounting for periodical boundary conditions) between the starting site of 
the monopole current at t = and the current position after i steps along the monopole trajectory, while L is the 
total length of the trajectory (all quantities measured in lattice spacing units). We have measured this quantity on 
trajectories wrapping only one time and tried to relate it to the monopole mass. We have considered the free particle 
approximation, thus definining a mass 

mFLUC = 2rL 2 "' (28) 

and obtaining the results reported in Table IVl and in Fig. El A better approximation would be to consider interactions 
at a mean field level. One can take the charge distribution around a monopole from the correlation function g(r) 
computed in Ref. Q and determine the potential around the position of a monopole induced by nearby monopoles and 
antimonopoles, then approximating it by an attractive harmonic potential: we have tried this procedure obtaining 
negligible corrections to the free particle approximation. 
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The mass obtained in this way shows some scaling deviations, especially at high T, which can be explained as 
follows. Wrapping monopole trajectories are not strict one-dimensional lines as in the one-particle path integral, 
but have some finite thickness, part of which is due to unphysical ultraviolet fluctuations in the form of monopole- 
antimonopole loops at the scale of the UV cutoff, which may be attached to a physical wrapping trajectory. Such 
unphysical fluctuations add to the physical ones that we are investigating. It is not easy to give a definition for a 
proper subtraction of the unphysical contribution. However we notice that it becomes negligible close to T c , where 
physical fluctuations of the trajectory become important: that is signalled by the fact that scaling is better and better 
as we approach T c . 

The value of the monopole mass obtained from spatial fluctuations of trajectories differs significantly from that 
obtained in the BEC approach: that is expected since tobec is related to the properties of multiple wrapping tra- 
jectories, while mpLUC to those of single wrapping trajectories: in both cases we have relied on some effective free 
particle behavior, but consistent results would have been obtained only if the particles were really free. It is however 
remarkable that also tofluc decreases rapidly as the transition temperature is approached, getting compatible with 
tobec if extrapolated down to T c (see Fig. |6|). 

Finally let us comment on Ax 2 measured on trajectories wrapping more than one time. In Table I VII we report 
the values of T 2 Ax 2 {k) measured for k = 1,2,3 at various temperatures. In the case of free particles one would 
expect Ax 2 (k) ~ kAx 2 (k — 1), instead, as anticipated in Section Till C[ values increase faster with k because of 
monopole-monopole repulsions. In particular Ax 2 (k = 2) and Ax 2 (k = 3) are respectively about 5-6 times and 10 
times larger than Ax 2 (k = 1). 

V. DISCUSSION AND CONCLUSIONS 

Let us summarize the main results obtained in this work. We have considered the statistical distribution in the 
number of wrappings in the Euclidean time direction, for thermal monopole trajectories exposed after MAG Abelian 
projection in the deconfincd phase of SU(2) pure gauge theory. We have shown that such distribution is compatible 
with that of a Bose gas condensing exactly at the confinement/deconfincmcnt critical temperature T c of the pure gauge 
theory. This result represents a strong link between confinement, monopole condensation and thermal monopoles. 
We have also found that in the approach to the criticality in the monopole ensemble, as T decrease to T c from above, 
the monopole masses are significantly reduced, and their interactions weaken. That fits into the overall idea of plasma 
gradually switching from electric to magnetic one. Our quantitative results about cluster densities shed further light 
on properties and interactions of the objects we call monopoles. Our investigation leaves some aspects that should be 
better understood and clarified and that we partially discuss in the following. 

1) We still do not understand the exact physical nature of thermal monopoles. The definition of monopoles studied 
in this work is strictly linked to a particular choice of Abelian projection, namely to that performed in Maximal 
Abelian Gauge, so that one should understand why some other Abelian projections give different, unphysical results 
and what non-Abelian gauge configurations are associated with the MAG Abelian monopoles. Some steps in this 
direction has been done in Ref. [TT| , where it has been shown a clear correlation of the non-Abelian gauge action density 
with the locations of thermal monopole trajectories. Objects identified by MAG projection have magnetic charge, yet 
various other quantum numbers are in principle possible. They can in principle be just monopoles, or dyons, or (in 
simulations with several quark flavors) they may also have various flavors, due to binding to quark- antiquark pairs 
or single quarks, or be a complicated mixture of all of them. Examples of fermionic excitations possessing magnetic 
charge and spin 1/2 are well known, and in fact required in any supersymmetric theory, complementing all monopole 
multiplcts. Scibcrg-Wittcn M = 2 superconformal theories with matter have not only magnetic monopoles with 
certain flavor allocations, but they also have intricate structure of dualities, relating several different (mutually dual) 
Lagrangians with the same theory. So, in simulations with quarks, one would need to separate Bose and Fermi-like 
behavior of various magnetic objects, which is by no means trivial. Those studies, as well as other issues related with 
the quest for understanding of the precise spectroscopy of the magnetic sector, will hopefully be continued elsewhere. 

2) While the statistical distribution of wrappings strongly suggests the presence of a BEC-like phenomenon asso- 
ciated with confinement, a clear description of the physical properties of the system has not yet been achieved. The 
ensemble of thermal monopoles is strongly interacting and its properties resemble those of a liquid [1, , therefore one 
should not try to fit its properties with those of a weakly interacting but rather with those of a strongly interacting 
Bose gas, like 4 He [241] . The effect of interactions is clearly visible in the typical shape of "polymers" made up of 
monopoles belonging to the same A;-cycle: we have not studied that in detail yet, but just remarked that from the 
scaling of the chemical potential we predict smooth linear shapes rather than the random walk behavior expected in 
the case of free particles; further investigation on this side could clarify more universal properties of the system. 

3) Another peculiarity of the system is that in fact it is a plasma of monopoles and antimonopolcs with Coulomb- 
like screened interactions, i.e. it resembles a Bose-Coulomb gas system which is neutral on average: in this sense the 
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chemical potential that we have introduced in our analysis should not be regarded as a true chemical potential, but 
rather as a parameter entering an effective description of the increase in the number of /c-cycles. 

4) The vanishing of the chemical potential ft, at the transition, with a critical exponent v 1 (see Eq. (|24|) ) compatible 
with the 3d Ising critical exponent v of the correlation length, shows that fi scales with the same critical behaviour of 
a mass. It would be interesting to better understand the meaning of this in the future, for instance by investigating 
the relation between the distribution in the number of wrappings of monopole trajectories and other quantities of dual 
nature, showing a critical behaviour above T c , like the dual string tension related to the spatial 't Hooft loop (33l - [35| . 

5) The fact that we still do not have a consistent physical description of the system manifests itself, for instance, 
in the absence of a consistent definition of effective monopole mass. However it is interesting to notice that different 
definitions tend to a common value m ~ T c when the condensation temperature is approached, where they are also 
roughly in agreement with the estimates reported in Ref. [24| . A similar physical scale comes out when studying 
the electric- magnetic asymmetry in Landau gauge (36[: it would be important to understand if the two scales are 
connected in some way. 

6) The fact that we obtain a condensation temperature Tbec = T c , within our statistical uncertainties, by just 
extrapolating data for the densities of trajectories with multiple wrappings, may be peculiar to SU(2) pure gauge 
theory, where the confinement/deconfinement transition is of second order nature, hence there is a critical behavior 
at the transition. In different cases, like for a first order transition, hence e.g. for the SU(3) pure gauge theory, while 
one expects monopole condensation to happen at T c anyway, the statistical analysis of the distribution of wrapping 
trajectories above T c could point to a condensation temperature slightly below the actual deconfincment temperature, 
i.e. to a sort of spinodal point. In this sense SU(2) pure gauge theory is a system where the connection between 
confinement and monopole condensation can be better established. We would like to clarify this point in the future 
by extending our investigation to SU(3) color gauge group, with and without dynamical quark contributions. 
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